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We investigate chiral and deconfinement transitions in the strong coupling lattice QCD for color 
SU(3). We combine the leading order Polyakov loop effective action of the strong coupling 
expansion and the next-to-next-to-leading order ( 1 /g 4 ) fermionic effective action with one species 
of unrooted staggered fermion. Two approximation schemes are adopted to evaluate the Polyakov 
loop effects; a Haar measure method (no fluctuation from the mean field) and a Weiss mean-field 
method (with fluctuations). The Polyakov loop is found to suppress the chiral condensate and 
to reduce the chiral transition temperature at jj. = 0. The chiral transition temperature roughly 
reproduces the Monte Carlo results in the region j3 = 2N c /g 2 < 4. 
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1. Introduction 

The phase transition in Quantum Chromodynamics (QCD) at finite temperature (T) and/or 
quark chemical potential (p) is attracting much attention in recent years. Since lattice Monte Carlo 
simulations have the notorious sign problem at finite p [Q], we need to invoke some approximations 
in QCD or effective models at finite chemical potential. Strong coupling lattice QCD (SC-LQCD) 
is one of the most promising approximation schemes in QCD to investigate the chiral and decon- 
finement phase transitions at finite T and p. The QCD phase transition has two characteristic 
features; the restoration of the chiral symmetry which is spontaneously broken in vacuum (chiral 
transition), and the liberation of color degrees of freedom which is confined at low temperatures 
(deconfinement transition). The chiral transition at finite T and p has been investigated in SC- 
LQCD [§, [||], ||]. The deconfinement transition at finite T was qualitatively explained on the basis 
of the mean-field treatment of the Polyakov loop in the leading order of the strong coupling expan- 
sion Higher order corrections on the Polyakov loop action has been investigated recently [Q]. 
One of the interesting developments describing the QCD phase transition at finite T and p can be 
found in the works g], which include the leading order Polyakov loop effective action (finite 
coupling constant) and the strong coupling limit effective action for quarks (infinite coupling con- 
stant), and enables us to describe the chiral and deconfinement transition in a single framework. 
To make a step forward towards the true phase diagram, it is necessary to consider Polyakov loop 
effects and finite coupling effects with quarks. 

In this proceedings, we develop & Polyakov loop extended strong coupling lattice QCD (P -SC- 
LQCD) framework JTot l by combining the leading order Polyakov loop effective action and NNLO 
quark effective action. We derive an analytic expression of the effective potential in P-SC-LQCD 
at finite T and pi, and investigate the chiral and deconfinement phase transitions at p = 0. 



2. Effective Potential with Polyakov Loop effects 

In the lattice QCD, the partition function and action with one species of unrooted staggered 
fermion for color SU(N C ) in the Euclidean spacetime are given as, 

iT LQ cD = J 9\X,X,U V ] e- Su *» , S LQ cd =S f + S g , (2.1) 

l d r i _ 2 

Z x p=0 x S p 

where x(2)> m o> U PtX , an d Up denote the quark (antiquark) field, bare quark mass, link variable, and 
plaquette, respectively. In the staggered fermion, the spinor index is reduced to the staggered phase 
factor r\ p x = (tJq^, f]j, x ) = , (— 1 )' toH Kc '~ 1 ) where p is the quark chemical potential. Note that 
we utilize the lattice unit a = 1 and represents physical values as dimensionless values normalized 
by the lattice spacing a throughout this proceedings. 

To treat the chiral and deconfinement phase transitions simultaneously, we evaluate the leading 
order Polyakov loop effects in the pure Yang-Mills sector and NNLO effects in the fermionic sector. 
In a finite T treatment of SC-LQCD, we first derive an effective action by carrying out spatial link 
(Uj) integrals, and evaluate the temporal link (Uo) integral later to consider the thermal effects of 
quarks [S]. 
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We can obtain the leading order Polyakov loop effective action S s from the spatial link inte- 
gral of the N T plaquettes sequentially connected in the temporal direction as shown in Fig. []]. 

where N T and P x = tr c Yl T Uo(z,x)/N c represent the temporal lattice size and Polyakov loop, re- 
spectively. This effective action shows the nearest-neighbor interaction of the Polyakov loop. The 
factor l/g 2 N c arises from the spatial link integral. By comparison, we derive NLO and NNLO 
quark effective action by the cluster expansion in Ref. [Q]. From a naive counting of the strong 




integral 



Figure 1: Leading order of the Polyakov loop effects in the strong coupling expansion. The squares in the 
left and loops in the right represent the temporal plaquettes and the Polyakov loops, respectively. 



coupling expansion order, the leading order Polyakov loop action composed of the plaquettes is in 
the higher order /g 2Nt )), compared with the NNLO terms (&(l/g 4 )) stem from the quark sec- 
tor. To describe the chiral and deconfinement phase transitions, we set the starting point of chiral 
and deconfinement dynamics as the effective action of ff(l/g°) and &{\jg T ) in the quark and 
Polyakov loop sectors, respectively. The Polyakov and NNLO quark effective action corresponds 
to an extension in the quark sector from this starting point. 

From the effective action described above, we obtain an approximate QCD partition function 
and define the effective potential J£" eff as 



"ell = - 7777 lo g^LQCD ^ "TT^log 



— SNNLO-^ e ff 



=& q (*;H,T) + U g {£,l) + , (2.4) 

where L and d{= 3) are the spatial lattice size and the spatial dimension, respectively. Here 

(x) 

SNNLO,^q,U g , and J^ff represent the NNLO quark effective action, quark free energy, pure glu- 
onic potential, and effective potential including only the auxiliary fields <J>, respectively. In this 
work, we first reduce the effective action to the bilinear form for the quark fields by introducing 
several auxiliary fields, and assume the mean-field values for <J> later. We obtain J^" q and U g by eval- 
uating the Grassmann {%,%) and temporal link (t/ ) integrals, In this proceedings, we evaluate the 
temporal link integral in two kinds of methods, the Haar measure and Weiss mean-field methods. 

We shall now derive the effective potential with Polyakov loop effects in the Haar measure 
method (H-method). In the H-method, we replace the Polyakov loop with a mean-field value 
and take into account the Haar measure in the Polyakov loop potential instead of carrying out the 
temporal link integral. The contribution to the effective action is, 

~-2p p L d tt, (2.5) 
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where fi p = (l/g 2 N c ) Nr N 2 d, £ = (P x ) and £ = (P x ). We assume the mean fields £ and £ are constant 
and isotropic. The temporal link integral is represented by using the Haar measure in the Polyakov 
gauge [0], which is a static and diagonalized gauge for temporal link variables, 



d% 



dldl-21 



l-6££ + 4(£ 3 + P)-3(££) 



(2.6) 



where ^b(x) = T7t Uq(x, t). This Haar measure shows the Jacobian in the transformation from the 
temporal link variables (Uq) to the Polyakov loop (£,£). jF q and U g are given as, 



-N c E q - TlogR(E q -fi,N c £,N c £) -TlogR (E q + p.,N c £,N c £) -N c \ogZ x , 



R{x,L,L) = 1 +Le- X I T +Le- 2x l T + e^ x l T , 
U g = -2Tp p ££ - Tlog [l - 6££ + 4 (£ 3 + 1 3 ) - 3 



(2.7) 
(2.8) 
(2.9) 



where fi p = {\/ g 2 N c ) l / T N 2 d, and Z? 9 is the quark excitation energy. Here we have replaced the 
N T with l/T and omitted irrelevant constants. J£" q includes the vacuum, quark, and antiquark free 
energies and the contribution of the wave function renormalization factor Z x . The quark free energy 
includes one- and two-quark excitations (e~^ E ~^l T ,e~ 2 ( E ~^l T ). In the confined phase {£ ~ £~ 0), 
the one- and two- quark excitations are suppressed and only the color-singlet state contributions 
remain. The pure gluonic potential U g does not include the fluctuation of the Polyakov loop since 
we treat the Polyakov loop as the mean field without the temporal link integral. In the Polyakov- 
loop extended Nambu-Jona-Lasino model, this pure gluonic potential is incorporated to express the 
properties of the deconfinement phase transition [|Tl|]. 

Now we shall evaluate the Polyakov loop effects in the Weiss mean-field method (W-method). 
The W-method includes some part of the fluctuation effects of the Polyakov loop. We first bosonize 
the Polyakov loop action by using the Extended Hubbard-Stratonovich (EHS) transformation [||], 
which is a procedure to bosonize the product of different types of composites. Then, the Polyakov 



loop action in Eq. (2.3) is linearized as, 



^eff 



g 2 N t 



N 2 £2| 

xj>0 



P x £ - £P X ) ~ 2P p L d ££ - 2p p £ (P x £ + £P X ) 



(2.10) 



where fi p ,P x , and P x are defined before. £ and /represent the auxiliary fields for the Polyakov loop, 
(£ = (P x },£ = {P x )). In "~" of Eq. ( 2.10j ), we assume constant and isotropic values for auxiliary 
fields £ and £. We obtain J£q and U g by evaluating the Grassmann {%,%) and temporal link (Uo) 
integrals, 



Z P 



Tlog(Z P /Lo)-N c logZ 



x i 



J^ndet 



2cosh(£ 9 /r) + ^ e^ /T + ^ V A/:r ] exp \rjtr^ + f]tr ^ c 



U g =2Tp p ££-TlogL l 



o • 



(2.11) 
(2.12) 
(2.13) 



We have defined tj = 2fi p £/N c and fj = 2fi p £/N c . In Eq. (2.12), we can perform the temporal 
link integral and obtain the analytic expression. Zp and Lq are functions of the combination of 
the modified Bessel functions. In Ref. [FTO], we show the explicit expression of them. Compared 
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with the effective potential in the H-method where the quark contributions are represented as the 
vacuum, quark and antiquark parts , Eq. (2.7), the effective potential is more complicated. 

Note that the pure gluonic potential U g includes the dependence on ljl explicitly (i.e. the 
dependence on When the quark chemical potential is zero (/I = 0), the Polyakov loop for quarks 
and antiquarks are the same {£ = £). In comparison, when the quark chemical potential is finite 
(jj. ^ 0), the Polyakov loop for antiquarks is generally different from that for quarks (I ^ £) [|I2|]. 



3. Chiral and Deconfinement Phase Transitions 




Figure 2: Left panel: Chiral condensate and Polyakov loop in P-SC-LQCD (solid lines), and chiral con- 
densate in SC-LQCD without the Polyakov loop effects (dashed line) as functions of T at fi = in the 
W-method. Right panel: Temperature dependence of do/dT (dashed line) and dl/dT (solid line) in the 
W-method. 



2 * ? (3 C (T=1/N T , u=0)(MC) • ' * 
2 T C (SCL, n=0)(MC) 




w/ PL (Haar) 

1 1 1 



1 2 3 . 4 5 6 

P=2N c /g 2 



Figure 3: Comparison of chiral transition temperature between P-SC-LQCD (solid line) and SC-LQCD 
without the Polyakov loop (dashed line). We show the results of Weiss mean-filed method (bolid solid line) 
and Haar measure method (thin solid line) in P-SC-LQCD. We define the critical temperature as the peak 
of —do/dT. The triangles represent the results of the critical temperature (T L .^=0, open triangle) and the 
critical coupling (j3 c , filled ttiangles) obtained in Monte-Carlo simulations with one species of unrooted 
staggered fermion [lj] 
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In Fig. [2], we show the chiral condensate (a), Polyakov loop, da/dT and d£/dT as functions 
of T. The chiral and deconfinement phase transitions seems to take place at almost the same T. 

We find that Polyakov loop suppresses the chiral condensate and reduces the chiral transition 
temperature in P-SC-LQCD. In left panel of Fig. ||] and Fig. [| we compare the chiral condensate 
o and the chiral transition temperature T Xfi with and without Polyakov loop effects. We find that 
both a and T xc in P-SC-LQCD are smaller than those without Polyakov loop effects; the chiral 
condensate becomes smaller when the Polyakov loop takes a finite value, and the chiral symmetry 
restoration takes place at lower T by the Polyakov loop effects. In SC-LQCD without Polyakov 
loop effects, we find the contribution only from color-singlet states, then we implicitly assume that 
the quarks are confined at any T. In the W-method, we have one- and two-quark contribution as 
well when the Polyakov loop takes a finite value. Quark excitation generally breaks the chiral 
condensates, then it promotes the chiral symmetry to be restored at lower T. A similar behavior is 
found also in the H-method. While qualitative behaviors are the same in both of the methods, the 
W-method exhibits a little larger T c than the H-method. Physically, the temporal link integral in the 
Weiss mean-field method favors the color-singlet states and therefore suppress the quark excitation. 

In our previous work on NLO and NNLO SC-LQCD []|], the critical temperature is calculated 
to be larger than MC results, and NNLO effects on T Xfi are found to be small. In Fig. we show 
T x c in two treatments of P-SC-LQCD in comparison with the MC results. The MC results, espe- 
cially in the region j8 < 4, are roughly explained in P-SC-LQCD developed in this work. Namely, 
the Polyakov loop and and finite coupling (NLO and NNLO) effects dominantly contribute to the 
reduction of T xc . This observation implies that introducing the Polyakov loop, the deconfinement 
order parameter, is essential to explain the QCD phase transition temperature. 

4. Concluding Remarks 

In this proceedings, we have derived an analytic expression of the effective potential at finite 
temperature and chemical potential in the strong coupling lattice QCD with the Polaykov loop 
effects using one species of unrooted staggered fermion. The chiral and deconfinement transitions 
at jJ. = are discussed with emphasis on the Polyakov loop effects on these transitions. The NNLO 
quark effective action in the strong coupling expansion is combined with the leading order Polyakov 
loop action, /g T ), and we have evaluated the temporal link (Uq) integral in two methods. One 
is the Haar measure method (H-method), where we replace the Polyakov loop with its constant 
mean-field without the Uq integral. The deconfinement dynamics is taken into account via the Haar 
measure. Another is the Weiss mean-field method (W-method), where we bosonize the effective 
action and carry out the temporal link integral explicitly, then the fluctuation effects are included. 
We have found that one- and two-quark excitations are allowed in both methods, when the Polyakov 
loop takes a finite value at high T. 

We find that the Polyakov loop reduces the chiral transition temperature and the obtained 
transition temperatures roughly explain the MC results in the region /3 < 4. The chiral and de- 
confinement transitions are found to take place at similar temperatures. The W- and H-methods 
exhibit qualitatively the same results, while the Polyakov loop effects are found to be weaker in 
the W-method. This is because the temporal link integral in the W-method favors the color-singlet 
states and leads to the suppression of the quark excitation. 
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As future works, we should investigate the QCD phase diagram in finite T and /I on the basis 
of the effective potentials derived in this work. Results in this direction are discussed in part in 
Ref. H. 
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